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ABSTRACT 

Por  a  genuinely  nonlinear  hyperbolic  system  of  conservation  laws  with 
added  artificial  viscosity,  -  u^  »  fin)--"  e  u^>  we  prove  that  traveling  wave 
profiles  for  small  amplitude  extreme  shocks  (the  slowest  and  fastest)  are 
linearly  stable  to  perturbations  in  initial  data  chosen  from  certain  spaces 
with  weighted  norm;  i.e. ,  we  show  that  the  spectrum  of  the  linearized  equation 
lies  strictly  in  the  left  half  plane,  except  for  a  simple  eigenvalue  at  the 
origin  (due  to  phase  translations  of  the  profile).  The  weight  -es~>is  used 
in  components  transverse  to  the  profile,  where,  for  an  extreme  shock,  the 
linearized  equation  is  dominated  by  unidirectional  convection. 
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SIGNIFICANCE  AND  EXPLANATION 


Many  squat Iona  of  mathematical  physics  take  the  form  of  nonlinear  systems 
of  conservation  laws.  Mathematically,  this  type  of  system  is  difficult  to 
deal  with!  generally,  smooth  solutions  must  develop  discontinuities  (shocks) 
in  finite  timei  thus  weak  (i.e.,  nondifferentiable)  solutions  must  be  admitted 
at  the  price  of  losing  uniqueness t  given  initial  data  may  be  continued  forward 
in  several  ways  in  the  class  of  weak  solutions.  An  extra  criterion,  called 
the  entropy  criterion  (from  the  fact  that  solutions  in  physics  have  increasing 
entropy)  is  imposed  on  the  solution  in  hopes  of  recovering  uniqueness.  One 
method  of  justifying  the  entropy  criterion  is  to  show  that  if  parabolic 
dissipation  terms  (e.g.,  viscosity),  generally  neglected,  are  taken  into 
account,  and  if  a  solution  to  the  inviscid  system  is  obtained  as  a  limit  of 
solutions  as  the  viscosity  vanishes,  then  this  solution  must  satisfy  the 
entropy  condition. 

Are  all  "entropy  solutions11  the  limit  of  solutions  of  the  viscous 
system?  So  far,  this  question  has  been  answered  only  for  scalar  equations, 
for  smooth  solutions,  and  for  simple  weak  solutiona,  namely  steady  shock 
waves.  Thus  a  steady  shock  satisfying  the  entropy  condition  is  the  limit  of 
smooth  "shock  profiles",  which  are  travelling  wave  solutions  of  the  associated 
parabolic  system. 

Insofar  as  the  incorporation  of  viscosity  into  the  system  yields  a  more 
accurate  physical  model,  shock  profiles  become  interesting  in  their  own 
right.  More  generally,  travelling  waves  of  parabolic  systems  arise  in  other 
important  contexts,  such  as  nerve  impulses  and  fronts  for  reaction- diffusion 
systems.  This  paper  is  concerned  with  the  stability  of  shock  profiles  with 
respect  to  perturbation,  with  fixed  viscosity  of  a  simple  but  nonphysical 
type.  As  apposed  to  the  situation  for  a  single  equation,  stability  results 
for  traveling  waves  of  nonlinear  parabolic  systems  are  few  and  far  between. 
However,  D.  Sattinger  has  established  a  linearised  stability  criterion  (one 
based  on  spectral  analysis)  for  such  waves,  which  can  insure  nonlinear 
stability  up  to  phase  shift.  Our  result  is  to  verify  this  criterion  for  a 
certain  class  of  perturbations  for  the  slowest  and  fastest  families  of  shock 
profiles  of  small  amplitude.  Unfortunately,  the  class  of  perturbations  we  can 
allow  cannot  be  treated  by  Sattinger *s  methods,  so  the  nonlinear  stability 
question  remains  open. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


LINEARIZED  STABILITY  OF  EXTREME  SHOCK  PROFILES  FOR 
SYSTEMS  OF  CONSERVATION  LAWS  WITH  VISCOSITY* 


Robert  L*  Pego 

*1.  Introduction 

In  this  paper  we  demonstrate  the  stability,  in  a  linearized  sense,  of 
viscous  shock  profiles  for  small  amplitude  extreme  shock  waves  of  a  hyperbolic 
system  of  conservation  laws. 

u.  +  f(u)  -  0  ,  u  e  R®  (1.D 

The  shock  profiles  are  travelling  wave  solutions  of  an  associated  parabolic 
system  obtained  by  adding  "artificial  viscosity", 

ufc  +  f(tt)x  *  WU^  (1.2) 

For  u  >  0  fixed,  we  exhibit  a  class  of  perturbations,  determined  by  a 
weighted  norm,  with  respect  to  which  the  profiles  satisfy  a  linearized 
stability  criterion  put  forth  by  Sattinger  (1976). 

We  assume  that  the  system  (1,1)  is  strictly  hyperbolic,  so  that  the 
matrix  df(u)  has  m  distinct  real  eigenvalues 

A,(u)  <  ...  <  Mu) 

with  corresponding  right  and  left  eigenvectors  r^Cu) ,  <^(u)  for  k  ■  1 

to  m,  with  *  r j  "  *  k-shock  for  the  system  (1.1)  is  a  two-valued 

weak  solution 


u(x,t)  - 


*L 

“r 


X  <  St 
X  >  St 


(1.3) 
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whara  u^  and  auat  satisfy  tha  Rankine-Hugoniot  jump  conditions 

fCu^)  -  flu^  “  •*UL“UR1  -  0 
and  tha  strict  entropy  inaqualitias 

\<0  >  *  >  W 
Vi"’*’  >  •  *  ViV 

Hera  X^(u)  is  assumed  to  ba  genuinely  nonlinear,  that  is  Vl^v^lu)  *  0. 

He  normalize  rk  so  that  71^  ^(u)  =  1 .  With  uR  fixed,  a  one 

He, 


parameter  family  of  k-shocks  exists,  with  ■  u  ( e)  and 
s  •  l^(uR) +  V2  e  for  e  >  0  small,  satisfying 


uNo) 


.“k 

du 

de 


(0) 


W  • 


An  elegant  proof  of  this  fact  may  be  found  in  Conlon  (1980).  Also  see 
Lax(1957). 

A  viscous  shock  profile  for  a  given  k- shock  (1.3)  is  a  smooth  traveling 
wave  solution  +((x-st)/ii)  of  the  system  (1.2)  such  that  +(€)  ♦  u^  as 
€  ♦  ♦(  £)  ♦  u^  as  C  ♦  +•.  The  profile  ♦(  5)  is  a  solution  of  a 

system  of  ODEs, 

-  f{#)  -  f(uR)  -  s(f*uR)  » 

joining  the  two  rest  point  u^  and  u^  .  In  the  present  situation, 

Foy  (1964)  proved  that  shock  profiles  $  (£  1  e)  exist 

for  small  amplitude  k-shocks  (with  s  ■  A^Ur)  +  hs).  His  proof  gives 
an  asymptotic  description  of  the  profile  as  e  -*■  0,  showing  that  to 
second  order  In  c,  the  profile  matches  a  properly  scaled  hyperbolic 
tangent  profile  for  the  scalar  Burgers'  equation,  ut  +  uux  *  uxx, 
along  the  eigenvector  rk(uR). 

We  are  concerned  with  the  asymptotic  stability  of  the  solution 
♦((x-st)/y;  c)  to  perturbations  In  the  Initial  data  for  a  fixed  p  >  0. 


Introduce  the  traveling  coordinate  £  ■  (x-st)/y  and  scale  t  by  t'  * 
t/y.  Then  (l/»  becomes 

■  u^  -  (df(u)  -  s)u^  .  (1.4) 

The  profile  $U;e)  is  a  stationary  solution  of  this  equation.  But 
so  Is  ♦U+y;  c)  for  any  phase  shift  y.  The  best  sort  of  stability  one  can 
expect  In  this  situation  Is  so-called  orbital  asymptotic  stability, 
which  means  that  Initial  data  of  the  form  u (e,0)  -  0  (5)  +  uQ  (5) 
yields  decay  of  the  form 

u(£,t)  -  $U  +  y)  +  0  as  t  -*•  » 

for  some  phase  shift  y,  If  u0  Is  small. 

The  approach  we  follow  Is  to  examine  the  linearized  stability 
of  the  profile  0.  Writing  u(5,t)  •  0(C)  +  v(C,t),  the  perturbation  v 
satisfies  a  nonlinear  evolution  equation 

vt  *  X(v)  .  (1.5) 

Linearizing  this  equation  at  v  -  0  (talcing  a  formal  Gateaux  derivative) 
we  get 

*  £  'v 

»  v^  -  (df($)  -  s)v^  -  (d2f (4>)4>^)v  .  (1.6) 

Typically,  If  one  can  show  that  the  spectrum  of  £'  Is  contained  strictly 
In  the  left  half  plane,  so  solutions  of  the  linearized  equation  decay 
exponentially,  then  one  can  show  that  small  solutions  of  the  nonlinear 
equation  also  decay  exponentially. 

However,  because  the  phase  of  the  profile  0  can  be  shifted  an  arbitrary 
-amount,  the  derivative  ^  Is  a  null  function  for£',  so  0  Is  In  the 
spectrum  of  £'.  This  situation  Is  typical  for  traveling  waves  In 
general,  which  arise  In  many  contexts  (e.g. ,  nerve  Impulses,  traveling 


fronts  for  reaction-diffusion  systems).  In  a  broad  setting,  Sattlnger 
(1976)  showed  nevertheless  that  If  0  Is  a  simple  eigenvalue  of  £  ',  and 

If  the  rest  of  the  spectrum  of  X*  on  a  suitable  function  space  Is  con¬ 
tained  strictly  In  the  left  half  plane,  then  Indeed  one  obtains  orbital 
asymptotic  stability  for  the  traveling  wave  $  with  exponential  adjust¬ 
ment.  Actually,  he  required  In  addition  that  the  resolvent  (A-X')*1 
satisfy  a  certain  asymptotic  estimate.  In  our  Appendix  we  show  that 
this  estimate  Is  automatically  valid  for  the  systems  he  considers. 

An  Important  feature  of  Sattinger's  analysis  is  the  use  of  spatially 
weighted  norms  for  function  spaces  to  "push"  the  essential  spectrum  (the 
spectrum  aside  from  Isolated  points  of  finite  multiplicity)  of  £’  to 
the  left.  It  may  be  verified,  In  fact,  that  the  essential  spectrum  of 
our  £  *  from  (1.6)  on  unweighted  Lp  spaces  Includes  the  origin  (see  the 
treatment  of  the  essential  spectrum  In  Henry  (1981),  appendix  to  Chapter 
5).  Our  main  result  Is  to  exhibit,  for  e  sufficiently  small  and  k  •  1 
or  m,  a  weighted  space  on  which  £•  satisfies  Sattlnger's  linearized 
stability  criteria.  Unfortunately,  our  result  does  not  Immediately 
yield  nonlinear  stability  by  Sattlnger's  theorem,  for  our  weight  fails 
to  satisfy  two  of  Sattlnger's  hypotheses.  In  particular,  (a)  the  weight 
Is  not  a  scalar  function,  and  (b)  it  Is  not  bounded  below  In  some  com¬ 
ponents,  so  the  nonlinear  terms  In  (1.5)  may  fall  to  be  continuous  on 
the  weighted  space.  (Xn  a  supplement,  we  describe  the  spectrum  on 
unweighted  spaces.) 

Our  main  result  only  applies  to  the  extreme  shock  profiles,  for 
which  k  *  1  or  m.  Our  approach  Is  to  try  to  decouple,  to  second  order 
In  c,  the  part  of  the  linearized  equation  (1.6)  along  r^u^)  from  the 
components  along  r^u^),  j  f  k.  Because  of  Foy's  result,  the  part  along 
r^  consists  mainly  of  Burgers'  equation  linearized  about  Its  scalar 


profile  (1  -  tanh  %x).  For  that  equation  the  weight  (cosh  %x)  Is 


standard.  The  components  along  r^  for  j  f  k  are  dominated  as  e  -*■  0 
by  convective  terms  -(Aj  -  .  For  an  extreme  shock  these  terms 

convect  In  the  same  direction.  Then  these  components  may  all  be  weighted 
the  same  to  achieve  exponential  decay,  which  Improves  as  e  «*  0  and 
permits  the  elimination  of  nonvanishing  "cross- terms"  coupling  the  rfc 
and  i-j  components  together. 

Let  us  now  state  our  main  result  for  1-shock  profiles  (m-shock 
profiles  may  be  treated  by  space  reversal,  x  +  -x).  By  a  linear  change 
of  coordinates,  we  may  assume 


df ( Ur)  *  dlag  (Xi(uR),...,  X^(uR))  . 

We  Introduce  the  following  spaces  of  functions  with  weighted  norms: 
For  c  >  0,  1  <  p  <  •,  define 

(Lp)*  »  {u:  F  -Fm  |  ulcosh  cx  €  Lp,  uje‘cx  €  LP,  j  -  2 . m  | 


with  norm 
Null 


_  «  max{llulcosh  cxll  _  ,  llu^e"cxll  _  ,  j  ■  2,...,m(  . 

*c  Lp  Lp 


Define  (C,,)™  and  II  II  _  similarly,  where  C„  Is  the  space  of  bounded 

U  t  II 

uniformly  continuous  functions  on  R  under  the  sup  norm. 


Theorem  1.1:  Fix  uR  €  Rm  and  suppose  Aj(u)  Is  genuinely  nonlinear. 
Consider  the  linearized  equation  (1.6),  vt  «  X'v,  for  the  evolution  of 
perturbations  of  the  1 -shock  profile  $(£;e).  Fix  c,  0  <  c  <  H-  Then 
there  exists  6,  0  <  6  <  c(l-c)  such  that  If  c  Is  sufficiently  small. 


the  spectrum  of  t '  on  any  of  the  spaces  (LP)jscC  ,  1  <  p  <  or  (Cu)J^c 
consists  of  a  simple  Isolated  eigenvalue  at  the  origin  and  a  part  which 
lies  In  a  sector  strictly  contained  In  the  left  half  plane* 

sa(_?<e2sJ  ■  U  6  C  ( Re( A  +  <  -(cos  o)|a  +  ke2& |}  , 

where  a  depends  on  e,  0  <  a  <  n/2. 

12.  Scaling  and  Transformation 

The  first  steps  In  our  proof  will  be:  (a)  to  scale  the  variables 
with  e,  using  Foy's  result  on  the  asymptotic  shape  of  the  profile  4>(£;e) 
as  e  -*•  0,  so  that  the  part  of  (1.6)  corresponding  to  the  linearized 
Burgers*  equation  appears  Invariant  of  e  as  e  ■*  0;  (b)  to  transform 
the  scaled  equation  by  Introducing  appropriate  weight  functions.  Theorem 
1.1  will  thereby  be  reduced  to  an  equivalent  statement  (Theorem  2.1)  con¬ 
cerned  with  the  spectrum  of  the  scaled  and  transformed  operator  on  an 
unweighted  space. 

Without  loss  of  generality,  assume  that  uR  «  0  *  f(uR).  Let  us 
recall  Foy's  asymptotic  description  of  the  1-shock  profile:  For  e  >  0 
sufficiently  small, 

♦U;e)  *  H efteZl  e) 
where  ty(x;c)  has  the  form 

<p(x;e)  *  <pB(x)r!(0)  +  e;(x;e)  , 

O 

where  4>  (x)  «  1  -  tanh  **x  and  sup  |c(x;e)|  <  C  independent  of  e. 

x 

(This  estimate  may  be  read  out  of  Foy's  proof  with  a  little  care.) 


Consider  the  terms  in  Eq.  (1.6).  By  Taylor's  theorem. 


£ — 


where 


df(4>U))  *  df (0)  +  d2f(0)*(O  ♦  R(0*2tt) 
d2fU(0)dg  *  d2f(0)^  +  S(0+(C)+g(g)  . 

1 

R(€)  “  /  (1-t)  dsf  (t#(0)  dt 


5(0  -  /  d*f(t#(0)  dt  . 

0 

R(E)  and  S U)  are  bounded  In  5,  and  for  each  ?  represent  a  trlllnear  map 
taking  a  vector  triple  to  a  vector.  The  equality 

(R4»2)^  *  5<M>£ 

holds,  where  each  side  Is  a  matrix. 

Let  us  now  Introduce  the  following  scaled  variables  (which,  despite 
appearances,  hav#  nothing  to  do  with  the  original  ones): 

x  ■  h  e£  .  t»  *  %e2t  ,  and  teu  *  v  . 


't  *  TV  •  VE  "  T  ux  •  *n<l  v«  *  T  uxx  • 


so  (1.6)  becomes 


ut  «  x*u  «  uxx  -  ♦  d2f  *(x)  ♦  >*efi(x/«sc)^(x)2J  ux 

-ld2f  ^x(x)  +  *se§  7,*c)f(xv  x(x)l  u  .  (2.1) 


i 


Here  the  derivatives  df  and  d2f  are  evaluated  at  0  (*  uR).  Define  the 
following  matrix-valued  functions: 

R(x)  *  d2f  c(x)  +  >sR(x/%e)  tp(x)2 
S(x)  *  d2f  ;x(x)  +  Js5(x(Jse)  ^(x)t|>x(x)  . 

Note  that  Rx  *  S,  and  R(x;e),  S(x;e)  are  uniformly  bounded  independent 
of  c.  We  also  define  the  matrix  A  *  (aj)  by 

A  *  d2f(0)M0)  .  so  a]  -  ti*d2f(r1,rj)  . 

It  is  true  that  a}  *  VXi*rx(0)  *  l(See  Lax,  1957).  Now  (2.1)  may  be 
written  in  a  slightly  more  coherent  fora: 

r'u  *  uxx  -  F(df  -  s)/*s€  +  dB(x)A  +  eR(x)J  ux  -  [*x(x)A  +  eRx(x)J  u 

(2.2) 

Observe  that  df-s  *  diag  ( ( A^-s)/»ae,  i  s  l,...,m)  and  that  ( A j  -  s)/*se 
■  -1  by  definition.  Let  us  write  (2.2)  in  a  convenient  block  form. 
Write 

u  *  |  j  ,  where  Q  *  (u2,...,um) 

Relative  to  this  decomposition,  also  write 


where  A  is  a  diagonal  matrix  with  positive  entries  X2-s,. . . ,Am-s. 
In  block  form,  then,  Eq.  (2.2)  is  written 


where  the  matrix  entries  are  operators  defined  by 

L’u1  *  LBu»  -  c[*y]x  ,  H'Q  *  -[(*B a.1  +  eR*)o]x 

Mu1  =  -  [(<}>Ba*i +  eRl)u1Jx  ,  Cu  *  LTD  -  £($BA  +  eR)u]  x  , 

B  T 

and  the  important  operators  L  and  L  are  defined  by 
lV  =  ujx  +  (1  -  *B(x))uJ  -  ^(xju1 
LTu  =  uxx  -  (AAe)ux  . 

This  completes  the  appropriate  scaling  of  the  variables  as  e  -*•  0. 

D 

The  operator  L  is  the  operator  that  would  be  obtained  by  linearizing 
Burgers'  equation  ut+uux  *  uxx  about  the  traveling  wave  solution 

D 

Q (x-t)  in  traveling  coordinates.  It  is  well  known  (Peletier,  1971; 
Sattinger,  1976)  that  considered  on  a  weighted  function  space  with 

D 

weight  w(x)  *  cosh  *sx,  the  operator  L  satisfies  the  linearized  stability 
criteria  we  are  concerned  with.  (Briefly,  the  idea  is  that  consideration 

D 

of  L  on  the  weighted  space  is  equivalent  to  consideration  of  the  trans- 
B 

formed  operator  wL  w"1  on  an  unweighted  space.)  Roughly,  what  we  would 

like  to  show, using  arguments  from  perturbation  theory  for  linear  operators, 
B 

is  that  L  "dominates''  the  first  component  of  £'  and  the  other  components 
only  "contribute"  spectrum  strictly  in  the  left  half  plane. 

Two  observations  are  of  primary  importance:  First,  in  general  a*  /  0, 
so  the  first  component  of  £'  apparently  cannot  be  considered  a  "small" 
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perturbation  of  L  .  Second,  observe  that  each  component  of  the  diagonal 
operator  iJ  consists  of  just  the  Laplaclan  plus  a  translational  term. 

The  time  evolution  for  this  operator  yields.  In  each  component,  solu¬ 
tions  of  the  heat  equation  In  a  frame  moving  at  constant  velocity.  No 
weight  which  is  bounded  below  can  yield  exponential  decay  rates  for  such 
solutions  (which  are  necessary  if  the  spectrum  of  LT  is  to  lie  strictly 
in  the  left  half  plane),  for  they  decay  to  zero  only  algebraically  in 
sup  norm.  However,  in  each  component  the  solutions  are  being  transported 
in  the  same  direction  at  increasingly  high  velocity  as  e  -*•  0  (since  we 
are  considering  an  extreme  shock).  Thus  one  may  expect  that  if  we  use 
a  weight  ecx  decaying  in  the  direction  of  translation,  we  would  obtain 
exponential  decay  for  solutions  of  u^  *  LTu  at  a  rate  which  improves 
as  e  -*■  0.  This  should  mean  that  the  spectrum  of  LT  moves  further  to 
the  left  as  e  0,  so  that  for  X  fixed,  one  expects  the  resolvent 
(X-lV  to  decay  in  norm  as  e  0.  Then  the  resolvent  equation 
(X-£')u  ■  f  might  be  solved  as  if  it  were  diagonally  dominant,  regard¬ 
less  of  the  coupling  terms. 

The  discussion  above  completely  describes  our  procedure,  and  moti¬ 
vates  the  following  transformation:  Fix  c  >  0,  and  in  the  block  form 
used  above,  define  the  matrix  weight 


Introduce  the  new  dependent  variable  z  ■  Wcu.  Then  the  equation  ut  * 
*'u  is  transformed  into 


ecx 
e‘cxlec* 

where  wc(x)  *  cosh  cx. 

In  this  way,  consideration  of  the  spectrum  of  the  original  £•  on 
P  m 

the  space  (L  )^ec  is  equivalent  to  consideration  of  the  spectrum  of 
£’c  on  the  unweighted  space  (LP)Jj  ,  except  for  a  scaling  factor.  To 
verify  the  scaling,  consider  the  following  equivalent  formulations  of 
the  resolvent  equation: 

(X  -  x')z  -  f(x)  In  ({.P)J  . 

Set  z  *  W cu,  f  *  UQg,  and  multiply  by  W”1  to  obtain 
(A  -  x')u  «  g(x)  In  (LP)J  . 

Now  set  v(0  ■  Heu(Jje£).  To  verify  the  scaling,  compute  formally 

i'v(C)  -  vt  *  Tut»  *  TT*r«*(x)  . 

So,  multiplying  by  e*/8,  we  obtain 

(he2 A  -£’)v  ■  h(0  In  (LP)J€C  , 

since  Wc(JseC)h(5)  e  (LP)q  .  Therefore,  A  Is  In  the  resolvent  set  of  ££ 
on  (LP)q  If  and  only  if  te2A  Is  In  the  resolvent  set  of  X'  on  0»P)£eC. 
Hence  Theorem  1.1  Is  equivalent  to: 


Zt  *  £>  *  Wc£'Wc-‘z 


w  Llw_1 
cL  c 


■cxmw:1 
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Theorem  2.1:  Fix  c,  0  <  c  <  H-  Then  there  exists  Bi  0  <  B  <  c(l-c), 
such  that  if  e  is  sufficiently  small,  then  the  spectrum  of  £'c  on  (LP)q 
1  <  p  <  oo,  or  (c  )".  consists  of  a  simple  isolated  eigenvalue  at  the 
origin  and  a  part  which  lies  in  the  sector 

Sa(-B)  =  1  A  e  (  |  fle(A  +  B)  <  -(cos  a)|A  +  B|  \ 

strictly  contained  in  the  left  half  plane,  where  a  may  depend  on  e, 
with  0  <  a  <  tr/2. 

The  proof  of  this  theorem  has  three  main  parts: 

(1)  Resolvent  estimates  for  the  operators  w/w'1  and  e"cxLTecx. 

C  t 

(2)  Verification  that  the  first  order  operator  w_H*ecx  has  uniformly 
bounded  coefficients.  It  suffices  to  show  that 

sup  |i|>(x;e)e2cxj  <  C  Independent  of  e. 
x 

This  yields  the  strict  requirement  c  <  %. 

(3)  Solution  of  the  resolvent  equation  fori'  by  block  Gaussian  ellmi- 
nation. 

In  the  next  two  sections  we  treat  (1).  The  estimate  for  lT,  Improv¬ 
ing  as  e  -*•  0,  Is  the  key  to  the  success  of  (3).  Since  iJ  Is  diagonal, 
it  suffices  to  treat  a  scalar  operator  uxx-2yux,  as  y  •*>  •.  Our  discus- 

D 

sion  of  L  parallels  Sattlnger's  treatment  (1976)  of  the  general  scalar 
equation.  The  argument  for  (2)  In  §5  consists  of  observing  that  the 
exponential  decay  obtained  In  a  standard  proof  of  the  stable  manifold 
theorem  holds  independent  of  c  as  e  ■+  0.  (i|>(x;e)  is  a  solution  of  an 

ODE  lying  in  the  stable  manifold  of  the  point  0.) 


»3.  Weighted  Resolvent  Estimates  for  uxx  -  2yux  as  y  ■*  « 


We  consider  the  scalar  differential  operator 


l,u  *  uxx  -  2yu, 


and  define 


L..U  «  e'cxL0ecxu  *  uyy  -  (2y  -  2c)u  -  (2yc  -  c2)u 

t  AA  A 


Observe  that  Lyu  ■  uxx-yau  Is  formally  self-adjoint.  We  are  concerned 
with  the  situation  y  >  c  >  0,  c  fixed,  y  large.  We  seek  norm  estimates 
Improving  as  y  «*-  •,  for  the  operators  (X-Lc)“l  and  (d/dx)(A  - 1^)”1 
on  Lp  and  Cu,  uniform  for  X  In  a  fixed  sector  of  C. 

The  following  result  suffices  for  our  purposes.  We  define 

Sa(B)  »  {X  e  C  |  «e(X-B)  >  -(cos  a) |X  -  b| }  , 

the  sector  with  vertex  B  €  R  and  angle  opening  2(ir-a)  to  the  right. 

Proposition  3.1:  Fix  B  €  R  and  c  >  0.  For  any  6  >  0,  If  y  Is 
sufficiently  large,  there  exists  a  <  ir/2  positive,  depending  on  y,  so 


x£Sa(e)  c 

for  j  «  0  and  1.  The  norm  Is  the  operator  norm  on  L*1  (or  Cu). 

This  proposition  follows  from  estimates  on  the  Green's  function 

for  X-L  ,  which  Is  obtained  by  transforming  through  1  .  The  resolvent 
c  y 

equation  for  Lc  Is  (A-Lc)u  ■  f.  Now 

(x-L.)  -  e<“-c>x(x- L..)e<c-u)x  , 


'  u  «  e*  t 


We  therefore  obtain  the  solution  formula. 


(X-Lcr‘f(x)  -  £  jr,-Y|x-y|*(u-c)(*-y)  f(y)  dy  , 


where  y(A)  *  V  y2  +  X  ,  «ey  >  0.  Thus  we  define 

ys.»)  •  £  +  <v-c)s  , 

so  (A-Lc)_1f  «  K^(«,A)*f.  When  «ey  >  y-c,  Ky(»,A)  €  L1  and  It  may  be 
verified  that  the  solution  formula  yields  a  bounded  Inverse  for  A-Lc 
on  Lp,  1  <  p  <  ,  on  the  domain 

0(LC)  *  |f  €  lp  |  f •  Is  absolutely  continuous,  f  »,f "  e  Lp  f 
(similarly  for  Cu  ).  We  omit  the  details,  since  a  similar  verlflca- 
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tlon  Is  carried  out  In  the  next  section  for  the  operator  *cL  w’  .  We 
also  find  that 

^(X-Lc)-’f(«)  ■  ^  .t-V+U-e>t*-y)f(y)  ^ 

^  mm 


So  defining 


dK^ts.x) 


dy]  . 


^  (-Y  ♦  y  -  c)  If  s  >  0 

i-Mu-c)  e(Y+y-c)s  If  s  <  0  , 


we  have 


d  (X-Lj'*f  »  dK(-,X)*f  . 


3x 


Using  Young's  Inequality,  HK*fN  _  <  IKK  Ifl,  we  find 

Lp  L1  Lp 


II  (X  -  Lc)“1 «  < 


and 


lliu-ij-l  «-i-( 

“  c  2WI  V 


_L  ( - i - ♦ - 1 - ) 

2|y|  \dtey  -  (y-c)  «ey  ♦  (y-c)  / 


j  /|y-(w-c)|  |y  ♦  (w-c)  |  \ 


2|y|  ^  flCY  -  (w-c)  <*eY  ♦  (w-c)  ) 


Proposition  3.2:  Let  y  *  yy*  +  X  for  X  e  |.  Then  If  y  >  2c, 
a  <  ir/2  may  be  chosen  so  that 
|y  -  (w-c) | 

(1)  -  <  347c  ,  If  X  €  S  (-yc) 

«*Y  -  (w-c)  a 

(2)  |y-  (w-c)  |  >  c/4  ,  1fX€Sa(-yc) 

Iy+(w-c)|  %  # 

(3)  -  <  tan  >i(ir-o)  <  2  for  X  €  SJ-y2)  . 

ter*  (w-c)  ° 


These  estlmtes  yield,  for  y  >  c,  X  €  $a(-yc), 

rfiH- 


Proposition  3.1  follows,  for  |y*  +  X|  >  (0  +  yJ)  sin  a  If  X  €  Sa(0), 
so  the  bounds  above  tend  to  zero  uniformly  for  X  €  S0(B)  «s  y  -*>  • 
at  the.  rate  y'*1. 
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Figure  4.  Action  of  the  Map  A  ■♦yy2  +  A 


Proof  of  Proposition  3.2:  Figure  4  gives  a  "before  and  after"  view 
of  the  complex  plane  under  the  map  X  -*•  ^y2  ♦  x  .  Part  (3)  of  the  propo¬ 
sition  becomes  obvious,  and  parts  (1)  and  (2)  are  implied  by  the  esti¬ 
mates 

(a)  sec  e  <  3.417c  and  (b)  y  >  (u-c)  ♦  c/4  . 

Calculation  of  y 
Define 

y  *  min  *  X 

A6S  J-cy) 


Me  Introduce  some  notation:  Set  0  ■  y2  -cy,  z  *  -cot  a  +  1  slna  , 


k  •  cos  a  •  s  >  sin  a.  For  t  >  0,  X  •  -yc  +  tz  lies  on  the  boundary 
of  S  (-cy)  and  y  -Vv1**  ■  /B*  tz  .  Me  minimize  «ey(t)  ■  «e/0+tz 
with  respect  to  t.  Now  (*ty)*  •  H(B-Kt*n(t)),  where  n(t)  -  |B  +  tz|. 
One  calculates  n'(t)  *  (t-B*0/|B  + tz|  .  Then 


[(fleY)1]’  ■  2(«ey)(«eY)' 


0 


only  when  |0  +  tz|ic  ■  t-8ic,  or 

ic2(02  -  2fkt  ♦  t2)  ■  t2  -  26wt  ♦  02ic2  , 

or,  since  we  seek  t  i  0,  t(l-w2)  ■  2BtcO-«2).  Therefore,  t  ■  20k,  whence 
(y)2  ■  Bs2.  so 

y  *  sin  oVy*  -  cy 

Now  B  ■  (y-c)2  ♦  c(y-c).  If  y/c  >  2,  then  c/(y-c)  <  1,  and  we  have 
>  (y-c) (l  +  {#•})  >  (v-0  ♦§ 


Hence  we  may  choose  a  <  v/Z  depending  on  y  so  that 

(sin  a)(y  *c*§)>(h-c*|-)  . 

Thus  y  >  (y-c)  +  c/4. 

Estimate  of  sec  6 
Define  6  by 


Iy  -  (y-c)| 

sec  0  »  max  - 

>eS  (-cy)  «rr  -  (y-c) 


Claim:  If  a  ■  ir/2,  then  sec20  ■  p/c  ♦  1.  Define 
h(t)  «  |y(t)-(p-c)|  «  | /6+tl  -  (p-c)| 

g(t)  »  «ey(t)  -  (y-c)  , 

where  6  *  P2-cp  as  above.  He  proceed  to  maximize  h2(t)/g2(t)  with 
respect  to  t.  This  quantity  Is  a  maximum  when  g2(h2)'  ■  h2(g2)\  We 
calculate  (with  n(t)  ■  |6+t1|): 

h2(t)  ■  |8  +  t1J  -  2(p-c)fley  ♦  (p-c) 2 
g2(t)  «  %(6+n(t))  -  2(p-c)«ey  +  (p-c)2 
(«ey)2'  «  2HeY(«eY)'  «  >*)r(t) 

(h2)'  -  n '( t)  -  2{p-c)n'(t)/4«eY 
(g2)'  •  %n»(t)  -  2(p-c)n'(t)/*teY 

Thus  h2/g2  Is  critical  when 

g2(2«ey  -  (prc))  ■  h2(«ey  -  (p-c))  »  h2g  , 
or 

gfley  ■  h2  -  g2 

Thus 

(Hey)2  -  (p-c)«ey  ■  H(n(t)  -  8)  , 
or 

*i(8 ♦  n(t))  -  (p-cHley  ■  »i(n(t)  -  8)  , 

$0  8*  (p-c)Aey.  Therefore,  p  ■  <tey  at  the  critical  point  t.  Squaring 
we  have  p2  ■  **(8  +  n(t)),  n(t)  ■  2p2  -8.  Squaring  again,  B2  ♦  t2  * 

4p*  -  4p28  ♦  82.  so  t2  ■  4p2(p2  -  8)  *  4p2(pc).  Therefore,  the  only 


r-’v 


W.i 


■W 

rr 

k' 


I 


j 


positive  critical  point  for  h2/g2  Is  t  ■  2p/uc.  We  calculate  the  value 
of  h2/g2  at  this  point: 

n(2p«^ic)  *  V02  ♦  4y*c  ■  Vy2C(u- c)2  +  4yc]  •  y(y  +  c)  . 

Since  <teY(2y«$c)  *  y,  we  have 

hi  «  v(v*c)  -  2(y-c)y  +  (y-c)2 _ 

g2  J*(y(y-c)  +  y(y+c))  -  2(y-c)y  +  (y-c)2 

t"2yttyc 

.  ctj_t_cj  .  k+  i  . 
c2  c 

At  t  ■  0,  h2/g2  •  1,  and  for  large  t,  h2/g2  <  2.  Therefore,  the  maximum 
value  of  h2/g2  Is  y/c  ♦  1,  establishing  our  claim  for  a  ■  ir/2. 

To  complete  estimate  (a)  for  sec  6,  fix  y,  and  observe  that 
h2/g2(t,a),  properly  defined,  1$  continuous  In  both  variables  so  long 
as  ir/2  -  o  Is  so  small  that  g  >  c/4.  Therefore,  If  ir/2  -  a  Is  suffi¬ 
ciently  small  (depending  on  y),  then 

max  ~  (t,o)  <  9  y/c  , 

0<t<»  g2 

establishing  (1)  of  Proposition  3.2. 


14.  The  Resolvent  Set  for  lBu  ■  uxx*  (tanh  %x)ux*  (>isech2>ix)u 

B 

In  this  section  we  discuss  the  resolvent  of  L  on  the  weighted  spaces 
(lp)*  and  (Cj!  with  weight  w_(x)  »  cosh  cx.  Equivalently,  we  study  the 

CMC  v 

operator 
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on  the  spaces  Lp,  1  <  p  <  «,  and  Cu.  L  may  be  obtained  by  linearizing 
Burgers'  equation  ut  +  uu  *  u  about  the  traveling  wave  solution  $Ux) 

«  1  -  tanh  He  in  traveling  coordinates.  Our  present  analysis  may  be 
considered  an  extension  of  that  of  Sattinger  (1976),  whose  results 

D 

apply  to  the  formally  self-adjoint  operator  Ly  Our  result  is  that 
if  0  <  c  <  H»  the  spectrum  of  L®  on  the  spaces  Lp  and  Cu  consists  of 
a  simple  eigenvalue  at  the  origin,  a  discrete  set  of  eigenvalues  in 
the  Interval  (-H,0),  and  a  part  contained  in  the  parabolic  region 

J  A  e  c  |  «eA  <  -c(l-c)  -  (yr2c)2  J  • 

Because  of  the  next  section,  our  results  for  c  strictly  less  than 
1/2  are  required  in  the  last  section.  In  conjunction  with  Sattinger* s 
treatment  of  nonlinear  stability,  these  results  imply  that  the  traveling 
wave  solution  of  Burgers'  equation  is  stable  to  sufficiently  small  per¬ 
turbations  which  decay  exponentially  as  |x|  -*■  ®,  no  matter  what  the 
rate.  A  similar  result  might  be  obtained  for  other  types  of  scalar 
equations  that  Sattinger  considers.  For  scalar  conservation  laws  with 
viscosity  however,  a  more  general  result  may  be  found  in  Il'in  and 
Oleinik  (I960)). 

Our  proof  uses  the  transformation  to  the  self-adjoint  form  to 
carry  information  about  that  well-studied  operator  (its  Green's  function 

the  location  and  multiplicity  of  Its  eigenvalues)  over  to  analyze  the 

B  B 

non-self-adjoint  1“  .  Let  us  collect  the  facts  we  require  about  , 

following  the  development  of  Sattinger  for  the  most  part. 

A  short  calculation  shows  that 


l?u  *  uv*  -  H(1  -  2  sech2%x)u 


Legna  4.1:  Define  y(\)  •  A  +  %  ,  «ey  >  0.  If  y  t  0.  the  homogen- 
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eous  equation  (l£-  X)u  *  0  has  a  system  of  solutions  $_•  ♦+,  and 
^(x(X)  with  the  following  asymptotic  properties: 

♦±  ■  e^O  +o(l))*  •  e±YX(±Y  +  o(l))  as  x  ♦  • 

t>±  •  •***(}  *ol D).  *±'  -  9±yx[±y*oO))  as  x  -►  —  . 

These  functions  are  single-valued  analytic  functions  of  X  In  the  complex 

plane  cut  from  -«•  to  along  the  real  axis. 

Proof:  This  Is  a  standard  result  on  the  asymptotic  behavior  of 

solutions  to  linear  ODEs.  To  obtain  seek  a  solution  of  the  form 

■  z_e'YX.  Then  z_  must  satisfy  the  equation  i*  •  2yz^  ♦  pz_  *  0, 

where  p(x)  ■  Ssech2%x.  Assuming  that  z_'  -*•  0  and  z_  -*■  1  is  x  •»  »  , 

Integrate  to  get  m 

zj(x)  •  J  pz.(s)  ds 

x 

*.(*)  •  1  -  /  (- — ^ - )p*.(s)  ds  • 

x 

Choose  x0  so  that 

/  p  ds  <  >i|y|  . 
x« 


Then  a  bounded  solution  z_(x)  may  be  obtained  on  the  Interval  [x0,») 
through  successive  approximations.  This  solution  does  Indeed  satisfy 
z  ■*  1,  zj  0  as  x  «*,  yielding  the  desired  behavior  for  $  . 


To  obtain  4>+,  set  $+  ■  z+e  .  Then  we  must  have  z+  +  2 yz+  +  pz+ 
*  0.  Set  z+(x0)  *  1,  z+(x0)  *  0.  Then,  integrating, 

z+(x)  *  “  /  e’^x”s^  pz+(s)  ds 
x0 


x0 


Our  choice  of  x0  ensures  that  this  last  equation  may  be  solved  for  a 
bounded  z+(x)  by  successive  approximations  on  the  Interval  [xo,00). 

This  solution  has  a  nonzero  limit  as  x  -►  “,  so  we  normalize  z+,  making 
tills  limit  1  and  obtain  $+  with  the  desired  asymptotic  properties. 

Because  of  symmetry,  we  may  set  ij>±(x)  ■  ^ (-x).  Since  z+  and 
z_(x,y)  are  obtained  by  successive  approximations,  they  are  analytic 
in  y  In  the  right  half  plane.  They  are  therefore  analytic  functions 
of  X  In  the  complex  plane  cut  from  -«  to  -h. 

Each  pair  or  clearly  forms  a  basis  of  solutions  to 
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the  homogeneous  equation  (L£  -X  )u  *  0.  Define  the  Wronsklan 

W(X)  ■ 

(It  Is  Independent  of  x).  The  homogeneous  equation  has  a  bounded  solu¬ 
tion  precisely  when  either  fley  ■  0  or  W(X)  ■  0.  If  «ey  *  0  (so  X  <  -fc), 
all  solutions  are  bounded.  If  W(x)  ■  0  but  Gey  >  0,  then  Is  a  mul¬ 
tiple  of  4>+,  which  then  spans  the  set  of  bounded  solutions.  Conversely, 

If  «ey  >  0  and  W(X)  f  0,  then  any  solution  bounded  as  x  -*■  »  Is  a  multiple 
of  but  ■  Aij»_  +  Bq/+,  where  A  f  0.  So  no  solution  can  remain  bounded 


for  all  x. 


Claim:  The  spectrum  of  L®  on  Lp  or  Cu  consists  exactly  of  the  zero 
set  of  W(x)  plus  the  Interval  (-»,-$*].  Clearly,  these  points  are  eigen- 
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values  for  L^  on  the  space  Cy.  If  W(A)  *  0  and  Gey  >  0,  then  the  eigen¬ 
function  ♦Jx.A)  Is  In  Lp  for  all  p  >  1.  If  Gey  *  0,  then  all  solutions 

of  the  homogeneous  equation  oscillate,  and  Is  bounded  but  not  In  Lp. 
It  Is  then  easy  to  verify  that  X  Is  In  the  approximate  point  spectrum 
of  L®  on  Lp,  1  <  p  <  ».  (Consider  the  sequence  un  ■  j(x/n)$_(x),  where 

j(x)  *  0  for  |x|  >1,  j (x)  *  1  for  |x|  <  «*,  and  j  Is  smooth.)  It  will 

o 

follow  from  our  analysis  of  the  resolvent  of  L®  below  that  X  Is  In  the 
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resolvent  set  for  when  Gey  >  0  and  W(A)  f  0. 

Lemma  4.2:  Let  X  be  a  zero  of  W(X)  In  the  complex  plane  cut  from 
-•  to  -h  along  the  real  axis.  Then  (a)  X  Is  real  and  Isolated,  and 
(b)  A  <  0. 

Proof:  (a)  $_(x,A)  decays  exponentially  as  x  +  +»,  so  defining 

00 

(u,v)  »  J  u(x)  v(x)  dx  , 

-00 

we  have 

X(<P_ )  ■  (lJ^,  , 

so  that  X  ■  X.  W(x)  Is  analytic  and  not  Identically  zero  In 
so  has  Isolated  zeros. 

(b)  The  function  $_(x,0)  *  Hsech  >$x  Is  positive.  (L®  kills  a 
positive  function  because  It  came  from  linearization  about  a  monotone 
traveling  wave.)  If  W(X)  *  0  for  some  X  >  0,  define 
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u(x)  -  4>.(x,X)/4>_(x,0)  . 

Then  since  y(X)  >  y(0),  u(x)  -*■  0  as  x  -*■  ±».  Also,  u  satisfies  the 
equation 

/  $'(x,0)  \ 

uwtf  +  (tanh  Hx  +  2 -  )  uw  -  Xu  *  0  . 

xx  V  *.(x,0)  /  x 

By  the  maximum  principle,  it  follows  that  u  =  0,  for  u  cannot  have  a 
positive  maximum  (at  which  uxx  <  0,  ux  *  0,  u  >  0),  nor  a  negative 
minimum.  (This  argument  applies  to  both  real  and  imaginary  parts  if  u 
is  not  real.)  This  contradicts  the  fact  u  t  0.  So  W(X)  f  0  for  X  >  0. 

O 

We  turn  now  to  construct  a  resolvent  formula  for  L".  Assume  that 
«ey  >  0  and  W(X)  f  0.  The  Green's  function  for  X  -  L®  is 

!— ]—  ^.(x,X)  (y,X)  for  x  <  y 

W(X) 

— —  <t>  (x,X)  ^.(y,X)  for  x  >  y  . 

W(x)  -  + 

D 

The  resolvent  formula  for  Lf  is 

*s 

00 

(X-L®)_1f(x)  =  J  K^x.y.X)  f(y)  dy  . 

•00 


Now  formally, 

(x-l \y' 


This  yields  the  resolvent  formula 


j  Kc(x,y,x)  f(y)  dy 


where 


Kc(x,y,X)  • 


^ ♦♦<*> fifrgr •  ♦» Ulr^  for  x<* 

^  ♦»  ifr§  •  *♦<*>  !§ir&  for  x>* 


D 

The  essential  step  In  proving  that  X  Is  In  the  resolvent  set  for  L" 
on  Lp  Is  to  prove  that  the  resolvent  formula  above  yields  a  bounded 
operator  on  Lp.  Here  Is  a  sufficient  condition  that  an  Integral 
operator  be  bounded  on  Lp: 


Lemma  4.3:  Let  K(x,y)  be  measurable,  with  II  K(x»* )li  <  C} 

L* 

Independent  of  x,  and  II  K(«  »y)ll Li  <  C2  Independent  of  y.  Then  the 

"«P 

*  -  /  K(xof)  f(y)  dy 

Is  bounded  on  Lp,  1  <  p  <  •,  with  nor*  »t  dost  C,1  " 


Proof:  Using  Jensen's  Inequality  and  Fublnl's  theorem,  we  find 


r 

J  /  Mx,y)f(y)  dy  dx 


00  00  /  00  ^ 

/  /  |f(y)lp  |Mx,y)j  dy  f  | K(x,s) |  dsj  dx 


•00  »oo 


w  w 

<  C? "V  |f(y)|P/  |K(x,y)|  dx  dy  <  Cp_1C2llfllPf 
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I 


The  following  estimate  Is  therefore  our  main  tool  for  describing  the 

fi 

resolvent  of  L": 

Proposition  4.4:  Assume  0  <  c  <  %,  flevT+T*  >  H-c  and  W(X)  f  0. 

Then  (a)  II  K  (x,*,X)ii  <  C(X)  Independent  of  x,  and 

c  l* 

(b)  II  K.(»  ,y,X)ll  <  C ( X )  Independent  of  y. 
c  L1 

Proof:  (a)  Assume  x  >  0  (the  estimate  for  x  <  0  Is  entirely 
similar).  We  estimate  the  Integral 


f  |Kc(x,y,X)|  dy 
•» 

separately  on  the  three  Intervals  (-*.0],  [0,x)t  and  [x,»).  On  the 
half  line  y  <  0,  and  for  x  >  0,  we  have  the  estimates 

♦♦MfSKJI  <  c, <*).<«"-<’»-'»*  '  (3.1) 

I'M*)  cosh  S  1  *  C,U).-<^*(%-c»*  *  C2(J)  (3 .2) 


Then,  since  x  >  0,  these  estimates  Imply 


f  |Kc(x,y,X)|  dy  < 


Cjl£i  _ _ 1_ 

|W(X)|  «CY-(>S-c) 


Now  Is  a  linear  combination  of  and  $+.  So  for  y  >  0,  we  estimate 
K<y>IfHy|  *  c»(>)  «(4teT  +  (,|‘c))y  .  (3.3) 


Together  with  the  estimate  (3.2),  this  yields 


/ |*c(x,y,A)|  dy  < 


C»C,  1 

I  W(A)  |  «ey+(i|-c) 


Finally,  on  the  half  line  y  >  0,  for  x  >  0  we  estimate 

!♦»!§&■§  |  *=  c.(x)  .*<«*  -  <--'»> 


(3.4) 


M*)fsr§|  *  c.(x)  «(<leY  *  (l,‘c))x  .  (3.5) 


which  yields  our  last  estimate. 


/  IM».y.x)|  4y  <  -£»E4-. - 3 - 

x  c  |M(»)I  «*y-0i-c) 


(b)  He  do  the  estimate  for 


/  |Kc(x,y,X)|  dx 

similarly.  Assume  that  y  >  0.  Then  for  x  <  0,  we  estimate 

koogtel  «  M*).** •  (3-«) 


Together  with  (3.4),  we  get  (since  y  >  0) 


f  l*c(x,y,A)|  dx  < 


ChC6  # _ } _ 

(W(X)I  «ey  -  (is -  c) 


The  estimates  (3,4)  and  (3.5)  yield 


/  |Kc(*.y.A)|  dx  < 


.lils _ 3 _ 

1«(X)(  «rr  -  («5-c) 
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Finally,  (3.2)  and  (3.3)  imply  that 


f  |Kc(x,y,X)|  dx 

y 


<  ...CjC,  1 

|W(A)|  fley  +  (*s-  c) 


The  estimate  for  y  <  0  is  similar. 


Remark:  Suppose  that  0  <  c  <  *s.  Then  fle/X  +  fc  >  %-c  if  and  only 
if 

«eX  >  -c(l-c)  . 

We  include  the  simple  proof:  fte/T  >  a  >  0  when  2(«e/F)2  *  tie  2  +  |z| 

>  2a2.  Then 

(«e  z)2  +  (im  z)2  >  (2a2  -ft  z)2  , 
so  ( Im  z)2  >  4a“-4a2«e  z  ,  and 

«ez  >  a2  -(^)*  . 

Now  take  z*A  +  %,  a  a  %-c. 

This  essentially  completes  the  determination  of  the  resolvent  set 

D 

for  Lc  .  To  complete  the  formal  proof  that  X  is  in  the  resolvent  set 
when  fley  >  *s-c  and  W(A)  i  0,  we  need  to  verify  that  the  (bounded) 
integral  operator  with  kernel  K  (x,y,A)  actually  yields  the  inverse 

w 

D 

of  X-Lc  on  a  suitable  domain.  This  verification  Is  straightforward; 
we  include  it  for  completeness. 

The  domain  of  L®  on  is  the  set  of  C1  functions  u  with  absolutely 
continuous  first  derivative  such  that  u,  u',  and  u"  are  in  Lp.  (On  Cu» 
require  u',  u"  in  C  . )  Given  f  in  Lp,  define 


-----  Y-j 


6f(x)  ■  J  Kc(x,y,A)  f(y)  dy  . 

•m 

ft 

To  show  that  6  Is  a  right  Inverse  for  A-i®  •  we  must  show  that  u  ■  Gf 
Is  In  the  domain  of  l®  and  that  (A-L®)u  «  f.  He  may  write 

x  • 

W(A)  ^  u(x)  •  *4(x)  J  im  ^  f(y)  dy  ♦  *Jx) J  ^  ^  f(y)  dy  . 

-*  x 

Clearly  u(x)  Is  absolutely  continuous,  so  may  be  differentiated  almost 
everywhere.  We  calculate 

x  • 

W(A)(^uj(x)  «  «x) y  ^  to  4  J  5T1  dy 

-•  x 

a.e.  Now  It  is  clear  that  u’  1$  absolutely  continuous.  One  may  verify 
that  u»  Is  In  Lp  using  the  proof  of  Proposition  3.4,  for  and  have 
the  same  asymptotic  properties  as  ^  and 
Differentiating  again,  we  find 

x 

W(X)(^U)(X)  "  ~*{X)  i£f(x)  *  *+(X)  f  ♦-  ^f(y)  ** 


♦  4j(x) 


m 

j  ^  dy 


Since  4.  and  are  eigenfunctions  of  ,  we  get 

a  *i.  w, 

(X‘V  w*  u  ■  ^  f  a.e. 


Therefore  (A-L®)u  ■  f  a.e. ..hence  In  Lp.  This  Implies  that  u  Is  In 

C 

the  domain  of  l®.  So  6  Is  a  right  Inverse  for  A-L®  on  Lp. 

B 

To  show  G  Is  a  left  Inverse  for  A-L"  on  Its  domain.  It  suffices 

D 

to  show  that  A  -  L®  Is  one  to  one.  (For  then  If  u  Is  In  the  domain,  and 

C 

f  ■  (A-L®)u,  we  have  (A-L®)(Gf-u)  ■  0,  so  Gf  *  u.)  But  If  A-L® 

C  C  v 

D 

Is  not  one  to  one,  there  exists  v  In  Its  domain  with  (A-Lc)v  ■  0. 

Then  v”  Is  absolutely  continuous,  so  v  Is  In  C2  and  Lp.  Since  A  Is 
not  an  eigenvalue  and  fley  >  c,  one  may  show  that  v  must  be  zero. 

Remark:  In  the  discussion  above  we  needed  to  know  that  If  A  Is  In 

B  B  _ 

the  resolvent  set  for  L®,  then  (d/dx)(A-  l“)  Is  a  bounded  operator 
on  Lp.  This  Is  actually  quite  a  general  fact  for  the  second  order  systems 

considered  In  Sattlnger  (1976).  Also,  In  our  treatment  we  did  not  follow 
Sattlnger  In  pursuing  resolvent  estimates  asymptotically  for  large 
| A| .  Estimates  such  as  he  requires  (see  Sattlnger  (1976),  Lemma 
3.411))  are  automatically  valid  for  the  class  of  second  order  (matrix) 
differential  operators  he  considers.  We  defer  a  discussion  of  these 
Issues  to  the  Appendix. 

Only  one  fact  remains  to  be  proved  In  this  section: 

Proposition  4.5:  For  0  <  c  <  H»  the  eigenvalue  0  Is  a  simple 
eigenvalue  for  l®  on  Lp  (and  Cu). 

Proof:  We  must  show  that  the  associated  projection  operator  on  Lp, 


Is  one-dimensional.  (Here  r  Is  •  smooth  closed  curve  In  the  region 

Oey  >  *i- c  enclosing  the  origin.) 

First,  observe  that  the  kernel  of  L"  Is  one-dimensional.  (If 
B  B 

L"u  ■  0,  where  u  Is  In  the  domain  of  L®,  then  u"  Is  absolutely  continu¬ 
ous,  so  u  Is  In  C*  and  Lp,  hence  Is  a  multiple  of  $_.)  He  claim  that 
the  quaslnllpotent  0£  associated  with  the  eigenvalue  zero, 

Dc  *  L5|,C  ■  EFT 1  d' 

r 

Is  actually  zero.  This  suffices  to  show  that  P£  has  one  dimensional 
range. 

We  argue  as  follows:  First  consider  the  case  c  «  %  on  the  Hilbert 

D 

space  L2.  The  operator  Is  then  actually  self-adjoint,  and  one 
obtains  »  0  on  L1  (we  cite  Kato  (1976),  V.3.5-6).  But  then  we 
may  write 

00 

Dcf(x)  *  2trT  /c  /  KC(x,y*^  dy  dc 

r 

00 

*  £  <»)  HT  A  /  Kis(x»y»;)  £  dy  dc  *  0 
r  -*  c 

for  any  f  In  C?(R).  But  Dc  Is  bounded  on  Lp,  so  Dc  »  0.  Therefore, 
the  eigenvalue  0  Is  simple. 


1 5.  Welghted_Est1mates  of  the  Shock  Profile 

The  k- shock  profile  $U»e)  obtained  from  Foy  (1964)  when  \{u)  Is 
genuinely  nonlinear  Is  :  edition  of  the  equation 

*  f(^)  *  s$  . 

(Here  s  *  Xk(uR)  +  He  and,  for  convenience,  uR  ■  0  «  f(uR).)  For  each 
e,  $U;e)  Is  In  the  stable  ran* fold  at  0  for  this  eqiiation.  We  need 
the  following  estimate,  for  the  1-snock  profile  in  particular): (see  the 
end  of  §2) 

Proposition:  Fix  c  <  H  positive.  Then  there  exists  a  constant  K, 
such  that  if  e  Is  sufficiently  small, 

sup  |<J>U;e)|  <  eKe”Ce€  . 

5 

Proof:  From  Foy's  proof,  4>(£;e)  =  Hei|>(HeE;e),  where 

i|>(x;e)  *  (1  -  tanh  Hx)rk(0)  +  et(x;e) 
and 

sup  |c(x;c) |  <  C  . 
x  €  F 
0  <  e<eo 

We  will  show  that 

sup  Mx;e)|  <  Ke"2cx 
x 

(independent  of  e).  ^  satisfies  the  equation 

*  Ae*  +  Fety)  (5.1) 

where 


and 


32 


Ae  •  diag  { (AjtO)  -  s)/He,  1  ■  . 

Then  Akk  *  -1  ,nd  0  *  f€<°>  "  dFC<°)-  Also*  <*2pe(v)  •  d2f(ev). 

Our  proof  follows  a  standard  proof  of  the  existence  of  the  stable 
manifold  based  on  solving  an  Integral  equation  by  successive  approxima¬ 
tions  (Coddlngton  and  Levinson  (1955)*  p.  330). 


In  block  form,  write 


where  A1  Is  diagonal  with  negative  eigenvalues.  A2  Is  diagonal  with  posi¬ 
tive  eigenvalues.  Choose  c  <  H  positive.  Then  If  e  Is  small  enough, 
there  Is  a  o  >  0  so  that 

|Ui(t) |  <  e^2c  +  0^T  for  t>0 
|U2(t)|  <  e0*  for  t  <  0 

Fixing  a  €  F*.  we  consider  the  Integral  equation 

T 

e(t,a)  •  Ui(x-T)a  ♦  f  lMt-s)Fc(e(s,a))  ds 

T 

m 

-  /  Ua(x-s)  Fc(6(s,a))  ds  .  (5.2) 

T 

The  following  Llpschltz  estimate  holds  for  Fc:  If  |6i|,|e2|  <  6,  then 

|Fe{e,)-Fc(ea)|  <  sup  |dF(e)|  |ea  - ea|  <  M6|e,-e2|  , 

|e|<6 

where  H  Is  a  bound  for  d2f  In  a  neighborhood  of  0. 
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Therefore,  restrict  6  and  a  so  that  2(a|  <  6  and  Hd(1/o  +  l/(2c+o)) 
<  We  solve  Eq.  (5.2)  by  successive  approximations:  set  0°(x,a)  =  0. 
Then  0!(x,a)  ■  Ui(t-T)a,  and 

|e*(T.«)-e*|  <  e‘2c(T'T)|a|  . 

By  Induction,  we  show  that 

|ei+1(T,a)  -  e4|  <  |a|  e‘2c(T*T)-2'A  . 

Indeed, 

|eU1(T,a) -0*|  <  /  e“^2c+0^T‘s^(M6|a|e"2c^s‘T^2"£+1 )  ds 

T 

♦  /  eo(T's)  (M6|a|e“2c(s-T).2-4+1)ds 

X 

-  M«|»|2-U,(e-2c(T-T)e-°T/  eos  ds 
'  T 

♦  eOT*2cT  j  ft-(2c+o)s  d  \ 

T 

<  H4|.|.-2^)  (I  +  ^)2-^ 

<  |a|e“2c^T‘T^2"£  . 

Therefore  0£  converges  uniformly  to  0(x,a)  with 

|e(T,a)|  <  2|a|e’2c(T"T)  for  x  >  T  .  (5.3) 


This  0  Is  a  solution  of  (5.1)  and  0(T,a)  has  a  special  form:  O^T.a) 


®j  If  j  <  k.  What  Is  more*  e(T,a)  is  the  unique  solution  of  (5.1) 
with  8j(T)  ■  8j  for  j  <  k  so  6  «*■  o  as  t  •*  *  (here  we  apply  the  stable 
manifold  theorem).  What  we  have  shown  Is  that  the  site  of  the  neigh* 
borhood  In  the  proof  of  the  stable  manifold  theorem  on  which  the  esti¬ 
mate  (5.3)  holds  does  not  shrink  to  zero  as  e  -►  0.  Instead,  this  estl 
mate  holds  whenever 


Now  i|»(r;e)  *  (1  -  tanh  T/2)rk  +  e;(T;e)»  where  |t|  <  C  Independent  of  e. 
Thus  If  eo  Is  sufficiently  small,  say  ea  <  5/4C,  then  there  exists  a 
fixed  T  so  that 


|i|»(T;e)|  <  6/2  for  all  e  <  e0 


Since  <i(t)  -*>  0  as  t  this  Implies  that  ^»(t)  *  0(t.i|>(T))  by  the 
stable  manifold  theorem  (for  e  fixed).  Therefore, 

Mt;c)|  <  (5eZcT)  e“2cr  for  all  e  <  e0  and  t  >  T  . 


Since  4»(x;e)  Is  bounded,  there  exists  K  such  that 

|^(t,c)|  <  Ke‘2cT  for  all  e  <  e0  and  t  real. 


§6.  The  Resolvent  Equation  for  l ^ 

We  proceed  to  complete  the  proof  of  Theorem  2.1  by  showing  that  If 
A  /  0  lies  exterior  to  a  suitable  sector  strictly  contained  In  the  left 
half  plane,  then  A  Is  in  the  resolvent  set  for  2"  ;  i.e.,  the  solution 
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operator  to  the  resolvent  equation  (jCc  -  X)u  »  f  is  a  bounded  operator 
on  (Lp)q  (or  fcu)-).  We  introduce  some  convenient  notation  for  the 
components  of  L  '  in  block  form  (see  §2): 


Also  set 


1  1  “1 
*  K  Lw 

Hi 

c  c 

c 

=  e'^w"1 

~Lc 

=  w  LBw~ 1 

R*z 

c  c 

c 

-  «-cx.  Tex 

-  e  L  e 

V2 

R‘r‘  -  -wJr'.w^V] 

-cxfHBr 


e  ~[(*BA  +  eR)ecxi]x  . 


Then 


K  ■ 


Lc  +  < 


L1  +  ft 
c  c 


The  operators  R*,  R  ,  and  N  are  first  order  differential  operators 

v  V  V 

with  smooth,  uniformly  bounded  coefficients  as  e  -*•  0.  Provided  c  <  h, 
the  same  is  true  for  because  of  the  result  of  the  previous  section: 


2  rv 

sup  |ip(x;e)  e  |  <  K  independent  of  e 
x 


Consider  the  resolvent  equation,  for  X  e  t  ,  f  €  (Lp)q  or  (cu)q: 


T,B 


Lc  +  eRc 


-  X 


Lc  +  *e 


-  X 


CHI) 


Fix  6  with  0  <  B  <  c(l-c)  and  a,  0  <  a  <  tt/2  so  that  the  only  point  in 


g 

Sa(-B)  (see  S3)  In  the  spectrum  of  L£  is  the  point  X  «  0.  Choose  r, 

0  <  r  <  6,  and  delete  from  the  right-facing  sector  Sa(-0)  the  disk  of 
radius  r  centered  at  the  origin,  obtaining  a  region 

pa  “  Sa(-B)  \  {  X  e  t  j  |X|  <  r  f  (see  Fig.  5)  . 

Claim:  If  e  is  sufficiently  small,  a  <  tt/2  may  be  chosen  (depending 

on  e)  so  that  P  is  contained  in  the  resolvent  set  of  £  '  . 

a  c 

Proof:  (1)  There  exists  C*  so  that  IIR*(L®  "  X)-I||<  Ci  for  all  X 
In  P  .  (For  this  we  require  some  asymptotic  control  on  ll(L"-X)  li  and 

U  L 

g  . 

II  (d/dx)(L°  -  X)"  ii  as  |X|  -*•  »  in  S^(-B) .  The  necessary  estimates,  proved 

w  CX 

in  detail  on  Cu,  are  found  In  the  Appendix.) 

(2)  For  e  <l/2Cj  ,  the  operator 

is  invertible  for  all  X  in  Pa,  with  norm 

H(L‘-xr1ll  <  2ll(L®-xr,||  . 

(3)  For  X  in  PQ  and  e  <  1/2C i  we  may  eliminate  N,.  from  the  resol¬ 
vent  equation,  obtaining 


(I  +  Q(lJ-X)_,)aJ-X) 


(')■(  "  ) 
V  z  /  Vf-Nj.aj-x)-^1  / 


where  Q  «  -  NC(L*  -  X)_1N*  . 

(4)  There  exists  C2  so  that  for  all  X  in  P 

1  a 
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There  exists 


IINCU>-  A)"1  II  <  C2  . 

(5)  Apply  Proposition  3.1  for  each  component  of  lI  . 

Co  >  0  so  that  if  c  <  e0t  we  may  find  o  <  n/2  so  that  if  X  e  p  ,  then 

I!NJ(lJ-  X)'1  ll  <  1/4C2 

and  T 

IIQ(L^-X)  *|  <  1/2  . 

Therefore,  we  may  Immediately  solve  for  z,  obtaining  the  estimate 

llzll  <  211  (L^  -  X)"1  II  (llfll  +  C2II  f Hi )  . 

Then  z1  Is  determined  by 

21  -  -  N*i) 

and  since 

i  -  (iJ-x)-'b(t -  sc(^-xr‘f‘)  , 

T  _  i  - 1 

where  B  *  (I  +  Q(LC-X)  )  ,  so  IIBII  <  2,  we  find  (using  (5)  to  estimate 

N*z)  that 

llz1!!  <  2II(L®-X)‘1||  (Ilf  lll  +  ^-  (llfll  +  Czllf1!! ))  . 

This  establishes  the  claim  above,  showing  the  Pa  is  in  the  resolvent 
set  of  2'c. 

To  complete  the  proof  of  Theorem  2.1,  it  remains  only  to  show  that 
inside  the  disk  |X|  <  r,  the  spectrum  of  £  '  consists  solely  of  a  simple 
eigenvalue  at  the  origin.  Consider  two  projections  P0  and  P£  defined 
on  (l*)*  or  (C„);  by 
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k:-- 


% 


In  block  form,  where 

PB  *  ST  /(5-l2>''  d« 

r 


Here  the  Integration  is  taken  in  the  positive  sense  around  a  circle  r 
centered  at  the  origin  with  radius  r,  r  <  ?  <  g.  From  the  theory  for 

D 

Lc  in  14,  P0  is  a  one-dimensional  projection. 


tr' 

£ 

K 

L-' 

F 

.  ** 

K-' 

C-; 

ft 


Claim:  If  e  Is  sufficiently  small,  then  l!Pe  -  P0 II  <  1. 

This  Implies  that  P  Is  one  dimensional.  But  £\  hence  £',  has 

c  C 

an  eigenvalue  at  0,  since  £  «  0.  The  point  zero  is  therefore  the 

only  point  Inside  the  circle  of  radius  r  in  the  spectrum  of  £ 

The  proof  of  this  claim  is  straightforward,  using  the  resolvent 
estimates  we  have  developed.  For  X  on  r,  f  in  (LP)J  ,  let 

( V )  ‘  (x>x)(f)  •  • 

Then  it  suffices  to  show  that  given  6  >  0,  then  for  e  sufficiently  small, 

llzll  <  fill f ll  and  Hz1  -  zBll  <  fill f II 

for  all  X  on  r.  The  first  inequality  is  clearly  guaranteed  by  our  pre- 

T 

vious  estimate  for  z  and  the  resolvent  estimates  for  L'  (Proposition 
3.1).  For  the  second,  consider  z1  in  more  detail: 
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Z1  »  (l®-a)‘,[i  +  cRc(l5-x)",J 

•  [fl  '  ^{lJ-X)“1B(?  -  Mc(^-xr,f1)]  . 

By  Proposition  3.1,  we  can  make  IIN*(lI- X)"1  li  as  small  as  we  like  by 

V  t 

choosing  e  sufficiently  small,  and  we  can  also  achieve 


II [i  ♦  cR*a®-x)_1]  1 


III  <  . 


from  which  we  may  extract  the  second  estimate.  So  the  claim  holds,  com¬ 
pleting  the  proof  of  Theorem  2.1. 


/ 

/7/ 


'K 


Figure  5.  The  Set  p 


APPENDIX: 

RESOLVENT  ESTIMATES  IN  WEIGHTED  NORM 
FOR  CERTAIN  PARABOLIC  SYSTEMS 


Here  we  develop  an  abstract  treatment  of  some  resolvent  estimates 
which  were  required  In  the  theory  of  asymptotic  orbital  stability  of 
traveling  waves  developed  by  Sattlnger  (1976).  For  the  class  of  opera¬ 
tors  and  spaces  with  scalar  weights  considered  by  Sattlnger,  we  show 
that  the  asymptotic  estimates  he  required  hold  automatically,  so  need  not  be 
separately  checked.  Similar  estimates  are  commonly  developed  for  para¬ 
bolic  operators  In  unweighted  spaces  In  the  theory  of  fractional  powers 
of  sectorial  operators,  which  are  generators  of  analytic  semigroups  (a 
good  reference  Is  Henry,  1981).  Our  treatment  Is  self-contained,  and 
proceeds  In  the  spirit  of  perturbation  theory  for  operators  generating 
quasi -bounded  semigroups.  At  the  end  of  this  section  we  prove  a  result 
Involving  matrix  weights,  validating  the  estimates  for  our  operator  £’ 
on  the  spaces  (CU)J!  considered  In  the  main  body  of  this  paper. 

Let  us  describe  the  estimates  Involved.  Consider  an  operator 

Lu  ■  Puxx  +  M(x)ux  +  N(x)u  , 

where  u(x)  €  Fm,  P  Is  a  positive  definite  matrix,  diagonal  for  simplicity, 
and  M(x)  and  N(x)  are  bounded  uniformly  continuous  matrix-valued  func¬ 
tions.  Let  w(x)  be  a  given  scalar  weight  function.  L  will  be  considered 
as  an  operator  on  the  weighted  space 

Bw  ■  )u:  F  Fm  |  w(x)u(x)  Is  bounded  and  uniformly  continuous  |  , 


I 

I 

I 

) 

I 

•v 

I 

* 


1 

r  - 


with  domain 

D(L)  ■  ju  €  Bw  |  ux  and  uxx  are  In  Bw  }  . 


The  space  Bw  Is  equipped  with  the  norm 

HullM  -  max  sup  |w(x)  u^x)! 
w  1  x 

We  also  define  a  norm 


nun 


W.l 


*  +  • 


The  weight  function  w(x)  should  satisfy 

(1)  w(x)  >  1  for  all  x; 

(11)  sup  |1  -  w(x+t)/w(x)|  ♦  0  as  t  0  . 
x 

In  Sattlnger's  framework,  L  was  obtained  by  linearizing  about  a  given 
traveling  wave  of  a  nonlinear  parabolic  system.  Condition  (1)  above 
Is  explicit,  and  (11)  Implicit,  In  Sattlnger's  analysis.  Condition 
(11)  ensures  that  the  shift  u(x)  ♦  u(x+t)  is  a  bounded  operator  on  Bw, 
continuous  In  t,  and  Implies  that  w(x)  Is  continuous  and  grows  only 
exponentially  as  |x|  -*■  ».  The  use  of  spaces  of  uniformly  continuous 
functions  Is  also  Implicit  In  Sattlnger's  work. 

The  second  proposition  below  contains  the  estimate  Sattlnger 
requires  In  Ms  Lemma  3.4  and  Theorem  4.1.  The  first  proposition  simply 
states  that  -l  Is  a  sectorial  operator  In  the  sense  of  Henry  (1981)  or 
m-sectorlal  In  the  sense  of  Kato  (1976). 


Proposition  A.l;  The  operator  L  on  Bw  with  domain  0(L)  Is  a  closed, 
densely  defined  operator.  For  some  a, 8  real,  with  0  <  a  <  ir/2,  the 
sector 
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{  A  €  t  I  «e(A  -  6)  >  -cos  a  |  A  -  B|  } 


Is  In  the  resolvent  set  of  L,  and  for  any  such  sector  S^, 


I  (A  -  L)-1f  II  w 


< 


l*-e| 


"f"w 


(A.l) 


for  all  A  €  Sog,  f  6  Bw. 

Proposition  A. 2:  For  any  A  In  the  resolvent  set  of  L,  ^  •  (A-  L) 
Is  a  bounded  operator  on  Bw,  and  for  any  sector  SQg  contained  In  the 
resolvent  set  of  L, 


MA-lpfl  Wf1 


< 


sS-j-  iif  ii 

a-bP 


(A.2) 


for  all  A  €  sag,  f  €  Bw. 


Our  approach  to  the  proofs  will  be  as  follows:  First  we  verify 
(A.l)  for  the  scalar  operator  u  -►  uxx;  then  (A.l)  Is  valid  for  the 
diagonal  operator  u  •+■  Puxx.  Me  then  establish  the  result  for  L  by 
treating  the  lower  order  terms  by  perturbation  arguments.  The  same 
procedure  Is  used  for  Proposition  A.2. 

Our  analysis  begins  with  a  study  of  the  translation  group  U(t) 
acting  on  Bw  by 

(U(t)u)(x)  «  u(x  +  t)  . 

Lemma  A. 3:  U(t)  Is  a  strongly  continuous,  quasi-bounded  group, 
meaning  that  for  some  constants  M  and  b, 

II U ( t ) u II  w  <  Heb't^llullw  (A. 3) 
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The  infinitesimal  generator  of  U(t)  is  the  operator 


Du  *  ux  , 

with  domain 

Z>(D)  =  { u  e  Bw  |  ux  e  Bw  f  . 

The  resolvent  set  of  D  includes  the  set  of  X  €  C  with  |<teX|  >  b,  and 
for  such  X  we  have  the  estimate 

II  (X  -  D)’1ull  <  - - -  llullw  .  (A. 4) 

w  |«eX|  -  b  " 

In  order  to  carry  out  the  proof,  we  first  estimate  the  weight: 

Claim:  There  exist  constants  M  and  b  such  that 

sup  w(x  +  t)/w(x)  <  Meb^  for  all  t  .  (A. 5) 

x 

Proof:  Using  property  (ii)  satisfied  by  w(x),  we  may  find  e  >  0 
and  B  >  0  so  that 

sup  sup  w(x  +  t)/w(x)  *  ebe  .  (A.6) 

1 1 1  <  e  x 

Given  any  real  t,  we  may  write  t  *  ne  +  l  ,  where  n  is  an  integer  and 
|t|  <  e.  Since  w(x  + je)/w(x+  (j  -  l)e)  <  ebe  for  any  j,  we  obtain 

w(x  +  t)/w(x)  <  (ebe)|n|«ebe  <  Meb|t|  , 

where  M  »  e2be.  So  the  claim  (A. 5)  is  established. 

Proof  of  Lemma  A. 3:  The  claim  above  implies  the  estimate  (A. 3). 
For  u  €  Bw,  let  us  verify  strong  continuity: 


sup  |w(x)(u{x  +  t)  -  u(x) ) | 

X 

<  sup  |(w(x)  -  w{x  + 1))  u(x  +  t)|  +  sup  |w(x  +  t)  u(x+t)  -  w(x)u(x)| 

X  X 

The  second  term  tends  to  zero  as  t  0  because  wu(x)  is  uniformly  con¬ 
tinuous,  and  the  first  term  tends  to  zero  because  wu(x  +  t)  is  bounded 
and 

sup  |1  -  w(x)/w(x  +  t)|  -*-0  as  t  *  0  . 
x 

Consider  the  operator  Du  *  ux.  It  is  not  hard  to  show  that  d[ D) 
is  dense  in  Bw>  We  claim  that,  if  flex  >  b,  then  X-D  Is  invertible, 
with 

CO 

(X  -  D)”lf  -  f  e"XtU(t)  f  dt  ,  f6Bw  .  (A. 7) 

0 

Indeed,  letting  v  denote  the  right-hand  side,  we  have 

o 

*v»w  <  /  e-t*tX.Mebt  dt  !«„  . 

0 

Also* 

v(x)  -  /  e‘Xtf(x  +  t)  dt  *  eXx  J  e_xt  f(t)  dt  , 

0  x 

so  v  is  differentiable  and  Xv-vx  ■  f,  so  v  Is  in  the  domain  of  D. 
Therefore  X-D  is  invertible  on  Bw  and  (A. 7)  holds. 

When  «eX  <  -b,  a  similar  analysis  holds  with  (A. 7)  replaced  by 


(X  -  D)-1f 


/  e’Xt  U(t)  f  dt 


(A.8) 


The  Identities  (A.7)  and  (A. 8)  Imply  that  D  Is  the  Infinitesimal  genera¬ 
tor  of  the  group  U(t).  but  we  will  not  use  this  fact.  Refer  to  Kato 
(1976)  for  the  details. 


Proof  of  Proposition  A.1  for  0*:  The  domain  of  the  operator  D2 
Is 

£>(D2)  «  Ju  e  Bw  I  ux  and  uxx  are  In  Bw  }  . 

Consider  the  resolvent  equation  for  D2: 

(X-D2)u  -  f  . 

Write  A  •  y2»  «ey  >  0.  Then  A-D2  •  (y-D)(y+D),  and  If  «ey  >  b,  then 
±Y  Is  In  the  resolvent  set  of  0.  so 


and 


u  -  (A  -  D2)-1f  -  (y  +  D)”1(y  •  D)”,f 


Hu«w  « 


- ^ -  Ifll 

(4teY-b)2 


So  the  resolvent  set  of  D2  Includes  the  region  * efi  >  b,  which  Is  "exterior 
to  the  parabola  fieA  «  b.  This  region  contains  a  sector  Sag.  where 
6  >  b2,  0  <  a  <  ir/2.  Let  Sog  now  denote  any  sector  In  the  resolvent 
set  of  D2,  with  0  <  o  <  v/2.  We  seek  to  show  (A.l)  holds.  The  resol¬ 
vent  Is  uniformly  bounded  In  any  compact  subset  of  SQg,  so  we  need  only 
verify  (A.l)  for  |A|  sufficiently  large.  It  suffices  to  prove  the  follow¬ 
ing: 


Claim:  There  exist  constants  C  and  c  so  that  If  A  €  Sae  and 
| A |  >  C,  then 


(«eA  -  b)2  >  c|X  -  B|  . 

For  the  proof,  first  choose  a,  0  <  5  <  a  and  pick  C  >  2|6|  so  large 
that  If  A  e  Sa6  and  JA|  >  C,  then  A  e  S~Q  ,  i.e.,  arg  A  <  it -5  .  Since 
the  set  «eA  ■  2b  is  a  parabola,  we  may  inflate  C  so  that  if  A  €  S  a 

ap 

and  | A |  >  C,  then  «e/\  >  2b.  We  also  have  arg  /k  <  (ir-a)/2,  so  it 
follows  that  («eA  -  b)/|  A|  has  a  positive  Infimum  ci  for  A  e  Sog, 

| A |  >  C.  Then 

(«evT  -  b)2  >  c 1 1 A | 

and  since  C  >  2|B|,  | A j  >  cs|A-b|  for  |A|  >  C,  c$  >  0,  and  the  claim 
follows.  Hence  Proposition  A.l  is  established  for  the  operator  D2. 

Now  Proposition  A.l  Is  valid  for  pD2,  p  >  0,  by  a  simple  scaling 
argument.  We  may  separately  analyze  each  component  of  the  resolvent 
equation  for  PD2, 

{A  -  PD2)u  *  f  , 

where  P  Is  a  positive  definite  diagonal  matrix,  and  find  that  Proposition 
A.l  Is  valid  for  this  operator  as  well. 

Our  treatment  of  the  lower  order  terms  of  L  is  based  on  the  notion 
of  relative  boundedness  of  closed  operators,  and  on  a  Landau- Kolmogorov 
inequality  for  generators  of  bounded  semigroups  (a  recent  reference  is 

Chemoff,  1979). 

Definition  (cf.  Kato,  1976):  Let  A  and  B  be  operators  on  a  Banach 
space,  with  D( B)  3  D( A).  B  is  said  to  be  A-bounded  with  bound  c0  if 
for  any  c  >  c0,  there  exists  C  so  that  for  all  v  €  p(A), 

HBvIl  <  Oil  vll  +  cIlAvll  . 
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Lemma  A. 4:  Let  M  and  b  be  given  from  Lemma  A. 3.  Then  for  any 
u  €  p(D2), 

It  ( D  -  b )  u  II  w  <  *ll(D-b)2ull  +  (M  + 1  )X  Hull  for  all  X>  0  (A. 9) 
ll(D-b)ull*  <  C0II  (D  -  b ) 2 u II  w  Hull  w  for  some  constant  C0  (A.  10) 

Corollary:  D  is  D2-bounded  with  bound  0  on  Bw. 

Proof:  For  X  >  0, 

II  (X  -  (D  -  b) )” 1 II  <  M/X  . 

For  u  €  d(D2),  we  may  write 

(X  +  D  -  b)u  «  (X- (D-b))-,(X2  -  (D-b)2)u  , 

so 

II  (D  -  b)ullw  <  XII  ullw  +  ”  (X2llullw  +  ll(D-b)2u|lw  , 

yielding  (A. 9).  Put  X  *  (II  (D-b)2ull ^/ll ull w)>5  to  obtain  (A.10).  For  the 
corollary  argue  as  follows:  0  is  (D-b)-bounded,  since 

II  Du  II w  <  ll(D-b)ullw  +  bllwl^  . 

By  (A. 9),  (D-b)  is  (D-b) 2 -bounded  with  bound  0.  Finally,  (D-b)2  is 
D2- bounded,  for 

II  (D- b ) 2 u II w  <  II D2 u II w  +  2bll  (D- b)uMw  +  b2ll ull w 

<  II D2 u II  +  2b(dl  (D  -  b ) 2 u II  ♦  C(e)llull)  +  b2llull  . 

™  WWW 

Taking  e  sufficiently  small,  we  find 


ll(D-b)2ullw  <  dl D2u II w  +  Cll  ull  w  .  (A. 11) 

Then  D  is  clearly  02 -bounded  with  bound  0. 

Proof  of  Proposition  A.l:  Let  Lu  s  Puxx  +  M(x)ux +  N(x)u,  and 
let  L0u  ■  Puxx,  Bu  *  M{x)ux+  N(x)u.  Since  H(x)  and  N(x)  are  bounded 
matrices,  8  is  D-bounded,  i.e. , 

H  Bu  It w  <  Cll  ux  II w  +  Qlullw  for  u  6  D(D)  . 

Applying  the  previous  corollary,  B  Is  D2 -bounded  with  bound  0,  hence 
Lo-bounded  with  bound  0.  But  then  L0  is  L-bounded,  for 


IILqU llw<  II Lull w  *  II Bull w  <  II Lull w  +  ellL0ullw  +  C(e)ll ull w  , 


and  If  e  <  1  we  obtain  II  L0u  II  w  <  cll  Lu  1^  +  Cll  ull  .  Hence  B  is  L-bounded 
with  bound  0. 

Fix  a  sector  Sag  with  0  <  a  <  n/2.  If  A  e  and  |A|  is  suffi¬ 
ciently  large,  A-L0  Is  Invertible  and  (A.l)  holds  for  L0.  Rewrite  (A.l) 
in  the  equivalent  form. 


"ull w  * 


-S*6—  ||  (A  -  L)u II 
I A  -  B| 


for  u  e  d( L)  . 


(A. 12) 


We  claim  that  (A. 12)  holds  for  A  in  Sag  if  |A|  is  sufficiently  large: 
For  u  €  d( L), 


II  (A  -  L  o )  u  II M  <  II  (A  -  L)ull  w  ♦  II  Bu  II  w  <  II  (A  -  L)u  II  w  +  cll  Lull  w  +  C(e)llull 
<  (1+e)  ll(A-L)ullw  +  (|A|e  +  C(e))  llullw  . 


Now  (A. 12)  holds  for  L0  for  some  constant  C  0.  Fix  e  <  1/4C  Q.  Then 


if  | X j  Is  sufficiently  large. 


— -2^-  (|X|e  +  C(e))  <  >s  , 
|X- 6| 


and  we  have 


l,u"w  « 


#  (X  -  L)u II  +  >1  lull 
|X-8|  w 


for  | X |  sufficiently  large,  yielding  (A. 12).  So  Proposition  A.l  is 
established. 


Proof  of  Proposition  A. 2:  Define  the  operators  L,  L0,  and  B  as 
above.  Let  us  show  that  If  A  Is  In  the  resolvent  set  of  L,  then 
D(X  -  L)"1  Is  a  bounded  operator  on  Bw.  Observe  that  since  D  Is  D2- 
bounded  with  bound  0,  and  L0  Is  L-bounded,  It  follows  that  D  Is  L- 
bounded  with  bound  0.  Then,  since  range  (A-L)”1  *  z?(L)  c  z?(D),  we 
have,  for  any  f  €  Bw, 

HD(X-L)‘,fllw  <  CllL(A-L)"lfllw  +  Cll  (X  -  L)‘lf  llw 

<  cll f II w  ♦  (c|Aj+C)  ■  (X -  L)“*f  II w 

<  C(A)  II f II w  . 


Let  us  now  assume  that  (A. 2)  holds  for  the  operator  L0  for  |X|  suffi¬ 
ciently  large,  A  In  any  given  sector  S0g,  0  <  a  <  tt/2,  and  show  that 
It  holds  for  L  If  |A|  Is  sufficiently  large.  Rewrite  (A.l)  In  the 
equivalent  form. 


lltflw  ♦  II Du II w  < 


H(A-L)ullw 


(A. 13] 


I 


for  u  e  d( L),  A  e  Sag,  | A J  sufficiently  large.  Now 
II  (A -  L o ) u II w  <  II  (A  -  L)u II w  +  HBullw 

<  ll(A-L)ullw  +  C(llull  +  II Dull )  . 

Since  (A. 13)  Is  assumed  to  hold  for  L0»  we  therefore  have 

(II  ull  w  +  II  Dull  w )  ^1  - 

Therefore,  If  |A|  Is  sufficiently  large,  (A.13)  holds  for  L  as  well, 
perhaps  with  a  larger  constant  C^. 

The  last  step  In  our  proof  Is  to  establish  (A.13)  for  the  operator 
D2.  By  considering  each  component  separately.  It  will  follow  Immediately 
that  (A.13)  Is  valid  for  the  diagonal  operator  L0.  Our  approach  Is 
use  the  Landau-Kolmogorov  Inequality  (A.10),  along  with  the  sectorial 
estimate  (A. 12)  for  D2. 

Fix  any  sector  Sag.  Then  for  A  In  this  sector,  |A|  sufficiently 
large,  A-D2  is  Invertible  and  (A.12)  holds.  Recall  from  (A. 11)  that 
(D-b)2  Is  D2-bounded.  Using  (A. 10),  we  find  that  for  u  e  p(D2), 

l(D-b)ulJ<  C0llullw(cllD2ullw  +  Cllullw)  . 

Now  ll  D2  u  II  w  <  II  (A  -  D2  )u  II  w  +  I A  |  II  ull  w  ,  so,  estimating  llullw  via  (A.12), 
we  obtain 

C  C  /  C  \ 

II  (0  -  b)u  II  *  <  — ^  (c  +  (c|A| +C) — SLi_\ I,  (A  -  D2)ull  *  . 

W  |  A  -  6|  \  |  A  -  0 1  /  w 

Therefore,  for  |A|  sufficiently  large. 


(A. 14) 


II  (D  -  b ) u  II  <  - SSi_  ||  (X  -  D2)u  II 

w  I*- el* 

Now  D  is  ( D-b ) -bounded,  and  we  have 


w 


II u II w  +  II Du II w  <  II  (D  -  b ) u II w  +  (b  +  l)llullw  . 

Using  (A. 14)  and  (A. 12),  we  obtain  (A. 13)  for  D2  for  all  US.  with 

dp 

1 A |  sufficiently  large.  This  concludes  the  proof  of  Proposition  A. 2. 


We  conclude  this  section  with  some  brief  remarks  about  matrix¬ 
valued  weight  functions  which  show  that  our  operator  t  •  of  (1.6)  satis¬ 
fies  the  estimates  (A.l)  and  (A. 2),  at  least  on  the  spaces  (Cu)^c 
for  e  sufficiently  small.  Consider  a  smooth  matrix  function  W(x)  which 
is  diagonal  for  all  x  for  simplicity,  with  W^x)  >  0  for  each  i. 

Define  the  weighted  space 

Bw  *  {  u:  F  -»•  Fm  |  Wu(x)  is  bounded  and  uniformly  continuous  }  , 

with  norm 

Hull  w  «  max  sup  fW.^x)  u^x)!  . 

For  example,  if  Wu(x)  *  cosh  cx  and  W..(x)  *  e"cx  for  j  >  1,  then 

J  J 

Bw  *  (Cu)*  .  The  theory  of  this  appendix  does  not  apply  in  By;  e.g., 
multiplication  by  a  constant,  nondiagonal  matrix  M  need  not  be  a  bounded 
operator  in  By.  However,  in  special  circumstances  the  estimates  (A.l) 
and  (A. 2)  may  hold  if  the  weight  W(x)  may  be  used  to  define  a  similarity 
transformation  which  takes  the  operator  under  consideration  to  a  "nice" 


-52- 


operator  acting  on  the  unweighted  space  (Cu) . 
Let  L  be  an  operator  on  Bu, 


Lu  *  Puxx  +  M(x)ux  +  N(x)u  , 

where  P  Is  a  positive  definite  diagonal  matrix,  and  M  and  N  are  smooth 
bounded  matrix-valued  functions.  Define  an  operator  on  (Cu)m  by 

LyV  ■  WLW” \  «  Pvxx  +  M(x)vx  +  N(x)v  . 

Now  M(x)  and  N(x),  In  general,  are  not  bounded  functions  of  x,  but  suppose 
they  are.  (This  Is  the  case  of  interest  in  this  paper,  where  Lw  corresponds 
to  ,  cf.  (2.4).)  Then  the  results  of  this  appendix  apply  for  Lw 
on  the  space  (Cu)m,  yielding  (A. 1 )  and  (A. 2)  for  Lw.  Consider  the 
equivalent  formu1at1ons(A.12)  and  (A.13).  Then  (A. 12)  holds  Immediately 
for  L  on  Bw*.  Given  any  u  €  D(l),  then  v  «  Wu  is  in  D(LW),  and 


II u  II u  »  llvll  <  — I  (X  -  WLW"1  )v  II  «  — 2S_ 

"  |A-e|  |  x  -  6| 


( X  -  L )  u  II  u  . 


To  verify  (A.13),  observe  that 

WDu  ■  Dv  -  Wxkf  *v 

Require  that  sup  | W  W” 1 (x) |  <  C  <  »  .  (This  Is  valid  for  the  weights  we 
x  * 

have  used  In  this  paper.)  Then 


II  u II u  +  II Du II w  <  (1  ♦C)llvHa#  +  II  Dv tl ^  < 


i»  -  si5 


(X-  Ly)vl 


l*-B| 


II  (X-  L)UllU  . 


Supplement:  Marginal  Stability  on  Unweighted  Spaces 

Here  we  show,  using  Theorem  1.1,  that  the  linearized  operator  £ '  for 
the  weak  1 -shock  profile  has  no  unstable  eigenvalues.  More  precisely  and  more 
generally: 

Theorem  For  e  >  0  sufficiently  small,  the  resolvent  set  of  £ '  on  the 
unweighted  space  (Cu)m  (also  (LP)B)  includes  all  complex  X  exterior  to  the 
sector  Sa(J/4t  3)  of  Theorem  1.1  and  exterior  to  some  parabolic  region  Pc  « 
{X|ReX  <  -C(Im  X)2}  which  lies  in  the  left  half  plane  but  includes  the 
origin. 

For  the  proof,  it  is  enough  to  show  that  £ '  has  no  eigenvalues  in  the 

region  described.  This  is  due  to  a  characterization  of  the  essential  spectrum 

(spectrum  aside  from  isolated  points  of  finite  multiplicity)  that  may  be  found 

in  Henry  (1981):  Let  A+  “  df(uR)  *•>*.“  df(u^)  ”  ■»  and  forra  the  sets 

S±  -  {X  e  $|  det(-t2  -  iT  A^  -  X)  «  0  for  some  real  T  }. 

Then  S+  U  s_  is  contained  in  the  essential  spectrum,  but  the  connected 

component  of  <:\ (S+  U  s_)  which  is  unbounded  to  the  right  is  devoid  of 

essential  spectrum.  Note  that  here,  any  X  €  is  of  the  form 

2 

Re  X  ■  -t  ,  Im  X  ■  t  (X  (u_  _ )  -  a)  for  some  j  -  1,...m, 

j  R#L 

2 

so  s+  U  s_  is  contained  in  a  parabolic  region  {X  €  $|Re  X  <  -C(Im  X)  }  for 

some  C  >  0.  Now  suppose  X  €  $  \  (S  (-V4  c20)  UP  )  and  X  is  an 

o’  c 

eigenvalue,  so  (X'-X)  ?  ■  0  for  some  function  T(x)  in  Dom( X* ) ,  (Therefore 
¥  is  bounded  and  smooth.)  It  will  be  convenient  not  to  scale  the 
independent  variable,  so  for  0  <  c  <  V2  ,  we  let 


We  claim  that  Wc¥  is  bounded,  so  that  X  is  an  eigenvalue  of  £^  , 
contradicting  Theorem  1.1. 


We  apply  standard 


Step  1.  We  claim  |e  '  ¥  (x)|  <  const,  x  <  0. 

theorems  on  asymptotic  behavior  of  solutions  of  linear  ODEs  (Coddington  and 
Levinson  1955).  As  x  +  f  is  asymptotic  to  a  solution  of  uKX  -  A.  ux  - 
Xu  ■  0,  i.e., 

¥  (x)  -  eW-X{ 1+o(1)),  and  ¥'(x)  -  p_eW-x(1+o<1 ) ) 

where  det  (p2  -  p_A_-X)  -  0  and  Re  p_  >  0.  In  this  case,  p_  is  of  the  . 

2  1/_ 

form  2p  -  (X.(u  )-s)  +  ((X.(uT)-s)  +  X)'2  for  some  j.  For  X  "exterior" 
-  3  L  j  L 

2  Vi 

to  S_,  we  have  Re( (X^ (uL)-s)  +  X)'2>  X^(uL)-s,  so  Re  p_  >  X^u^-s.  Recall 

.  2 
that  s  -  X  (u  )  +  '2e'  *  1 '  and  ul  *  UR  +  er1  +  0(  e  >#  80 

•% 

Re  p  >  e  “  0(e).  For  e  small,  we  then  must  have  Re  p_  -  ec/2  >0, 
so  the  claim  holds. 

Step  2  We  claim  |eCcx/^2  f 1 |  <  const,  x  >  0.  A  different  argument  is  needed 

1  * 

here.  In  block  form,  ¥  •  (¥  ,  ¥),  and 

(L1-!) (cosh  ecx/2  f1]  -  -M1  [tTtcx/2f] 
c  c 


The  right  side  is  bounded,  from  step  1.  By  the  construction  of  the  resolvent 

of  £'  in  §  6,  X  is  in  the  resolvent  set  of  L1  if  X  ?  0  is  exterior  to 
c  c 

2  1 

S  (-e  8/4 ) .  Therefore  we  may  find  4  in  Dom(L  ),  so  bounded,  with 
Q  C 


Hence 


,  1  , .  .  1 .  - ecx/2 “ 

(L  -X)  ♦  -  -M  [e  ¥] . 

c  c 

(L^-X)^1  -sech  ecx/2  ♦]  ■  0 


The  expression  in  brackets  is  bounded,  so  has  asymptotic  behavior 
eM+x(1+o(1))  as  x  ♦  +•,  where  (since  X^(uR)-s  ■  -e/2) 
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